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■ Abstract Accurate quantum Monte Carlo calculations of ground states and low- 
lying excited states of light p-shell nuclei are now possible for realistic nuclear Hamil- 
tonians that fit nucleon-nucleon scattering data. Results for more than 30 different 
(J n ; T ) states, plus isobaric analogs, in A < 8 nuclei have been obtained with an excel- 
lent reproduction of the experimental energy spectrum. These microscopic calculations 
show that nuclear structure, including both single-particle and clustering aspects, can 
be explained starting from elementary two- and three-nucleon interactions. Various 
density and momentum distributions, electromagnetic form factors, and spectroscopic 
factors have also been computed, as well as electroweak capture reactions of astro- 
physical interest. 
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1. INTRODUCTION 


A major goal of nuclear physics is to understand the stability, structure, and re- 
actions of nuclei as a consequence of the interactions among individual nucleons. 
This extremely challenging many-body problem is exacerbated by the fact that we 


*The US government has the right to retain a nonexclusive, royalty-free license in and to 
any copyright covering this paper. 


53 



Annu. Rev. Nucl. Part. Sci. 2001.51:53-90. Downloaded from arjoumals.annualreviews.org 
by 193.219.53.38 on 09/26/05. For personal use only. 


54 


PIEPER ■ WIRINGA 


do not know in detail what the interactions are. Quantum chromodynamics, the 
fundamental theory for strong interactions, is so difficult to solve in the nonper- 
turbative regime of low-energy nuclear physics that no one is currently able to 
quantitatively describe the interactions between two nucleons from first princi- 
ples. However, we do have a tremendous amount of experimental data on nucleon- 
nucleon (AW) scattering, and we can construct accurate representations of AW 
interactions by means of two-body potentials. These potentials are complicated, 
depending on the relative positions, spins, isospins, and momenta of the nucleons, 
which makes finding accurate quantum-mechanical solutions of many-nucleon 
bound and scattering states a demanding task. Nevertheless, a variety of precise 
calculations have been carried out for the three- and four-nucleon systems, which 
reveal an additional complication: Two-nucleon forces alone are inadequate to 
quantitatively explain either the bound-state properties or a variety of Nd and Nt 
scattering data. This is not surprising, because meson-exchange theory and the 
composite nature of nucleons indicate that many-nucleon forces, at least three- 
nucleon forces, are important. However, these are even more difficult to construct 
from first principles. In like manner, meson exchange introduces two-body charge 
and current operators, and calculations of electroweak transitions show that these 
are also necessary to achieve agreement with data. 

Despite these difficulties, tremendous progress has been made in the past 
decade, both in the characterization of nuclear forces and currents and in the de- 
velopment of accurate many-body techniques for evaluating them. A multienergy 
partial-wave analysis of elastic AW scattering data below T \ a b = 350 MeV, pro- 
duced by the Nijmegen group in 1993 (1), demonstrated that over 4300 data points 
could be fit with a/ 2 per datum ~ 1. This was achieved by careful attention 
to the known long-range part of the AW interaction and a rigorous winnowing of 
inconsistent data. Their work inspired the construction of a number of new AW 
potential models that gave comparably good fits to the data base. 

Concurrently, several essentially exact methods have been developed for the 
study of few-nucleon (A = 3, 4) systems with such realistic AW potentials. The 
first accurate, converged calculations of three-nucleon (NNN) bound states were 
obtained by configuration- space Faddeev methods in 1985 (2), and the first accu- 
rate 4 He bound states were obtained by the Green’s function Monte Carlo method 
in 1987 (3). Other methods of comparable accuracy have been developed since 
then for both A — 3, 4 bound states and NNN scattering. A comprehensive re- 
view of many aspects of the few-nucleon problem can be found in the article by 
Carlson & Schiavilla (4). The result of all the progress in the past decade is that 
we have obtained an excellent understanding of few-nucleon systems, including 
both structure and reactions, although a few unanswered questions remain. 

Until recently, nuclei larger than A = 4 were described within the framework 
of the nuclear shell model or of mean-field pictures such as Skyrme models. The 
classic description of p-shell (5 < A < 16) nuclei was developed by Cohen & 
Kurath in the 1960s (5). In that work, a large number of experimental levels were 
fit in terms of a relatively small number of two-body matrix elements, but no 
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direct connection was made to the underlying NN force. Beginning with the work 
of Kuo & Brown (6), approximate methods for relating the shell-model matrix 
elements to the underlying forces were developed and applied to even larger nuclei. 
However, in most calculations, only the valence nucleons, i.e., the nucleons in the 
last unfilled shell, are treated as active, whereas the core is inert. Recently, the 
first “no-core” shell-model calculations, in which all the nucleons are active, have 
been made in the p-shell with effective interactions that have been derived by a 
G-matrix procedure from a realistic underlying NN force (7). These calculations 
agree with the other exact methods for the A = 3, 4 nuclei, but fully converged 
calculations for larger systems have not yet been obtained. 

The present article focuses on recent developments in quantum Monte Carlo 
methods that are making the light p-shell nuclei accessible at a level of accuracy 
close to that obtained for the s-shell (A < 4) nuclei. The quantum Monte Carlo 
methods include variational Monte Carlo (VMC) and Green’s function Monte 
Carlo (GFMC) methods. VMC is an approximate method (8-1 1) that is used as 
a starting point for the more accurate GFMC calculations. Because of its relative 
simplicity, VMC has also been used to make first studies of various reactions, 
such as the nuclear response to electron (12, 13) or pion (14) scattering, and for 
electroweak capture reactions of astrophysical interest (15, 16). The GFMC method 
is exact in the sense that binding energies can be calculated with an accuracy of 
better than 2%; it has been used to compute ground states and low-lying excitations 
of all the A < 8 nuclei (3, 17—19). These studies of light p-shell nuclei have shown 
that nuclear structure can be predicted from the elementary nuclear forces. As 
presently implemented, both the VMC and GFMC methods require significant 
additional computer resources for every nucleon added, so they probably will be 
limited to systems of A < 12 nuclei. A new method, whose accuracy lies between 
that of the VMC and GFMC methods, is auxiliary-field diffusion Monte Carlo 
(AFDMC); it has been implemented for large pure neutron systems (20) and with 
further development should be capable of handling nuclei above A — 12. 

Most of the calculations discussed here have been made using Hamiltonians 
containing the Argonne tqg NN potential (21), alone or with one of the Urbana (9) 
or Illinois (22) series of NNN potentials. Argonne v i g (AVI 8) is representative of 
the modern NN potentials that give accurate fits to scattering data; the Urbana and 
Illinois models are modern NNN potentials based on meson-exchange interactions 
and fit to the binding energies of light nuclei. These models are described in some 
detail in Section 2 to show the complexity of modern nuclear forces and the conse- 
quent challenge for many -body theory. The VMC method is presented in Section 3; 
it starts with the construction of a variational wave function of specified angular 
momentum, parity, and isospin, 'Ty(7 7r ; T), using products of two- and three -body 
correlation operators acting on a fully antisymmetrized set of one-body basis states 
coupled to specific quantum numbers. Metropolis Monte Carlo integration (23) is 
used to evaluate ( 'T y I H | TV ) , giving an upper bound to the energy of the state . The 
GFMC method, described in Section 4, is a stochastic method that systematically 
improves on TV by projecting out excited-state contamination using the Euclidean 
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propagation 'T(r) = exp [— (H — E( l )r]'R l /. In the limit r -> oo, this leads to the 
exact ( H ). 

Energy results for the ground states and many low-lying excited states of the 
A < 8 nuclei are presented in Section 5. With the new Hamiltonians, the excitation 
spectra are reproduced very well, including the relative stability of different nu- 
clei and the splittings between excited states. Energy differences between isobaric 
analog states and isospin-mixing matrix elements have also been computed to 
study the charge-independence breaking of the strong and electromagnetic forces. 
A variety of density and momentum distributions, electromagnetic moments and 
form factors, transition densities, spectroscopic factors, and intrinsic deforma- 
tions are presented in Section 6. One of the chief benefits of this approach is 
that absolute normalizations can be computed, and effective charges so commonly 
needed in traditional shell-model calculations of electroweak transitions are not 
required. 

Although the quantum Monte Carlo methods have been used primarily for sta- 
tionary states, it is also possible to study scattering states, as discussed in Section 7. 
Some first calculations of astrophysical radiative capture reactions in the p-shell 
are presented in Section 8. Both of these applications are in their early stages but 
promise to be major fields of endeavor for future work. In Section 9, we discuss 
neutron drops, systems of interacting neutrons confined in an artificial external 
well, that can serve as reference points for the development of Skyrme interactions 
for use in very large neutron-rich nuclei. Finally, in Section 10, we consider the 
outlook for future developments. 


2. HAMILTONIAN AND CURRENTS 


The Hamiltonian used in this review includes nonrelativistic kinetic energy, the 
AV18 AW potential (21), and either the Urbana IX (17) or one of the new Illinois 
series of NNN potentials (22): 

H = £*'■ + £%+ E v &- 1. 

i i<j i<j<k 


The kinetic energy operator has charge-independent (Cl) and charge-symmetry- 
breaking (CSB) components, the latter due to the difference between proton and 
neutron masses. 


K t = K 


ci 


K? SB = — IVf - -I — - — ItrfVA 


2 . 


where t, is the isospin of nucleon i. 

The AVI 8 potential is written as a sum of electromagnetic and one-pion- 
exchange terms and a shorter-range phenomenological part. 


Vii = V: + V ", + vf : 


3. 
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The electromagnetic terms include one- and two-photon-exchange Coulomb in- 
teraction, vacuum polarization, Darwin-Foldy, and magnetic moment terms, all 
with appropriate proton and neutron form factors. The one-pion-exchange part of 
the potential includes the charge-dependent (CD) terms due to the difference in 
neutral and charged pion masses. It can be written as 

v U = f2 (^) \ f nc 2 [X l] r l ■ Tj + 4. 

where 7) , = 3 r zl r ZJ — r, ■ r ; is the isotensor operator and 

Xij=\(X? j + 2 X±), 5. 

= 6 - 

X'I’j — [Y (mr ij)cr i ■ (Tj + T(mr IJ )S i j]. 7. 

Here Y (x) and T(x) are the normal Yukawa and tensor functions: 

Y(x) — — £y(r), 8. 

x 

T(x) = (^+ j + 1 )n*)?r(r), 9. 

with short-range cutoff functions £y(r) and £-/ (r), and the X Jl are calculated with 
m = m n ± and m n o. The coupling constant in Equation 4 is f 2 — 0.075 and the 
scaling mass m s — m n ±. 

The remaining terms are of intermediate (two-pion-exchange) and short range, 
with some 40 adjustable parameters. The one-pion-exchange and the remaining 
phenomenological part of the potential can be written as a sum of 18 operators, 

V U + V U= J2 v P^ij)°ij- 10 - 

p=l,18 

The first 14 are Cl operators. 

Off 1 ' 14 - [1, cr, • (Tj. Sij, L • S, L 2 , L 2 (<t, • a-), (L • S) 2 ] ® [1, r, • r ; ], 11. 

while the last four, 

Oj’j =l5A& = [1, (Ti ■ (Tj, Sij] ® Tjj, and (r zi + r zj ), 12. 

are CD and CSB terms. The potential was fitted directly to the Nijmegen NN 
scattering data base (1) containing 1787 pp and 2514 np data in the range 0-350 
MeV, to the nn scattering length, and to the deuteron binding energy, with a/ 2 per 
datum of 1.09. 

Although the CD and CSB terms are small, there is a clear need for their 
presence. The Nijmegen group studied many AW potentials from the 1980s and 
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Figure 1 Three-body-force Feynman diagrams: a is the Fujita-Miyazawa or two- 
pion P-wave, b is the two-pion S-wave, and c and d are three-pion rings with one-A 
intermediate states. 


before, and found that potentials fitted to np data in 7’ = 1 states did not fit pp data 
well even after allowing for standard Coulomb effects, and vice versa (24). All 
five modern AW potentials, Reid93, Nijm I and II (25), CD Bonn (26), and AV18, 
which fit pp and np scattering data with a / 2 /point near 1, have CD components. 
The CSB term is required to fit the difference between the pp and nn scattering 
lengths and is consistent with the mass difference between 3 H and 3 He, as shown 
below. The identification of the proper Cl, CD, and CSB AW force components is 
important in setting the proper baseline for the introduction of NNN forces, which 
are required to fit the 3 H and 3 He binding energies as a primary constraint. 

The Urbana series of three-nucleon potentials is written as a sum of two-pion- 
exchange P-wave and remaining shorter-range phenomenological terms. 


Vijk - v% p + Vj* . 


13. 


The structure of the two-pion P-wave exchange term with an intermediate A exci- 
tation (Figure la) was originally written down by Fujita & Miyazawa (27); it can 
be expressed simply as 

Vm P = E KrfXy- X *k){Ti ■ Tj • Tj ■ r,} 

eye 

+ C 2 P n [Xf j ,X^][T i -T j ,rj-r k ]), 14. 


where Xf- is constructed with an average pion mass, m n = \m n o + | m„±, and 
^T cyc is a sum over the three cyclic exchanges of nucleons i, j, k. For the Urbana 
models — j A , as in the original Fujita-Miyazawa model (27), whereas 
other potentials, such as the Tucson-Melbourne (28) and Brazil (29) models, have 
a ratio slightly larger than The shorter-range phenomenological term is given 
by 


V; 


ijk 


= E a 


RT 2 (m n rij)T 2 (m n r jk ). 


15. 


eye 


For the Urbana IX (UIX) model, the parameters were determined by fitting the 
binding energy of 3 H and the density of nuclear matter in conjunction with AVI 8. 
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As shown below, although the combined AV18/UIX Hamiltonian reproduces 
the binding energies of s-shell nuclei, it does not do so well for light p-shell nuclei. 
Recently a new class of NNN potentials, called the Illinois models, has been devel- 
oped to address this problem (22). Five models, IL1-IL5, have been introduced; 
in this review, we concentrate on IL2. These potentials contain the Urbana terms 
and two additional terms, resulting in a total of four coupling constants that can 
be adjusted to fit the data. 

One term, is due to ixN S-wave scattering as illustrated in Figure lb. It 

has been included in a number of NNN potentials such as the Tucson-Melbourne 
(28) and Brazil (29) models. The Illinois models use the form recommended in the 
latest Texas model (30), where chiral symmetry is used to constrain the structure 
of the interaction. However, in practice, this term is much smaller than the V?%’ P 
contribution and behaves similarly in light nuclei, so it is difficult to establish its 
strength independently just from calculations of energy levels. 

A more important addition is a simplified form for three-pion rings containing 
one or two A resonances (Figure lc, d). As discussed elsewhere (22), these dia- 
grams result in a large number of terms, the most important of which are used to 
construct the Illinois models: 


A R 

V ijk 


- A ar 

— A 3 7T 


50 , , 26 , 

-s'X + jKa 


i 

<7 


16. 


Here S* and ,4 1 are operators that are symmetric or antisymmetric under any 
exchange of the three nucleons, and the subscript x = a or r indicates that the 
operators act on, respectively, spin and space or just isospin degrees of freedom. 

The S' T is a projector onto isospin- 1 triples: 

, 2 

S T — 2 + -(T; -Tj +Tj ■ T k + T k ' Ti) = 4P r=3/2 . 17. 

To the extent isospin is conserved, there are no such triples in the s-shell nuclei, 
and so this term does not affect them. It is also zero for Nd scattering. However, the 
S ' T term is attractive in all the p-shell nuclei studied. The A l T is the same isospin 
commutator, [t, ■ Tj. Tj ■ Tk], that appears in V^ ,P (Equation 14), but the A^A^ 
term is repulsive in all nuclei studied so far. In p-shell nuclei, the magnitude of the 
A^A^ term is smaller than that of the term, so the net effect of the V^’ AR is 

slight repulsion in s-shell nuclei and larger attraction in p-shell nuclei. The reader 
is referred to Reference 22 and its appendix for the complete structure of V ( ^’ AS . 

Relativistic corrections to the Hamiltonian of Equation 1 have been studied in 
three- and four-body nuclei (31-33). We only briefly outline the results here. If a 
square-root kinetic energy. 


K* el = 


p] + mj - m, 


18. 


where p, is the momentum of nucleon i, is introduced into the Hamiltonian, then 
the AWpotential must be readjusted to refit the AWdata. This can be done for AVI 8 
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with only small adjustments of the parameters, to produce a u rd . VMC calculations 
show that the resulting (K rel + v rel > ~ (K + v) for 3 H and 4 He; thus the square-root 
kinetic energy may be neglected in light nuclei. 

A second relativistic correction, the so-called boost correction, arises from the 
fact that the NN interaction u l; - depends both on the relative momentum p = 
(p, — Pjj/2 and the total momentum P = p, + p ; of the interacting nucleons. 


Vij( p, P) = %(P) + Svijip, P), 19. 

where Svjj( p, P = 0) = 0. Fitting the two-body potential to NN scattering data 
determines only t>, ; . The Svjj can be determined from D, ; - ; a suitable approximation 
(31-35) is 

• r p - v _ »,j], 20 . 

where only the first six terms (the static terms) of Vjj are retained and gradients 
of tensor operators are ignored. The 8v,j is meant to be used only in first-order 
perturbation in a nonrelativistic wave function. Its expectation values in light nuclei 
are nearly proportional to (VA,.) (22), and thus one can consider that part of ,4 R , 
Equation 15, comes from Sv/j. This proportionality does not, however, extend to 
nuclear matter or pure neutron systems. There, one must use a reduced A* R that 
does not contain the part ascribable to 8 Vjj , and explicitly evaluate {8 i;, ; ) ; doing 
so results in a softer equation of state for dense matter. 


3. VARIATIONAL MONTE CARLO 

The variational method can be used to obtain approximate solutions to the many- 
body Schrodinger equation, // V = FT, for a wide range of nuclear systems, 
including few-body nuclei, light closed-shell nuclei, nuclear matter, and neutron 
stars (36). A suitably parameterized wave function, TV, is used to calculate an 
upper bound to the exact ground-state energy. 


(Ty|//|Ty) 

(S'vl'M 


> Eq. 


21 . 


The parameters in TV are varied to minimize Ey, and the lowest value is taken 
as the approximate solution. Upper bounds to excited states are also obtainable if 
they have different quantum numbers from the ground state, or from small-basis 
diagonalizations if they have the same quantum numbers. The corresponding TV 
can be used to calculate other properties, such as particle density or electromagnetic 
moments, or to start a Green’s function Monte Carlo calculation. In this section we 
describe the Ansatz for TV for light nuclei and briefly review how the expectation 
value is evaluated and the parameters of TV are fixed. 
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The best variational wave function has the form (11) 


I'M 



E 0 u ljk + ufjn + Eu, ff 

i<j<k i<j 


I*P). 


where the pair wave function, T P , is given by 


22 . 


|4M = 






23. 


The Ujj , U(j S , Uijk - and U™ 1 are noncommuting NN and NNN correlation oper- 
ators, and S is a symmetrization operator. The form of the antisymmetric Jastrow 
wave function, T j, depends on the nuclear state under investigation. For the s-shell 
nuclei, the simple form 


|W,> = 


n fuk( r ik,r jk )Ylfc(nj) 


i< j<k 


l< t> A (JMTT})) 


24. 


is used. Here / e (r ly ) and f^ k are central (spin-isospin independent) two- and three- 
body correlation functions and <t> A is an antisymmetrized spin-isospin state, e.g., 

® 3 (nn)) = ■ 4lt p 4 p f n)l 25 ' 

|<f>4(0000)) =A\'\p\p / \n] f n), 26. 


with A the antisymmetrization operator. 

The two-body correlation operators 6', ; and !//■ s are sums of spin, isospin, 
tensor, and spin-orbit terms: 


U, 


= E 

P= 2,6 L k^i,j 


n fuk (r ’> 


, r ik ) 


u„(r,;)0' 


27. 


U. 


LS 


= E 

P= 7,8 L kjtij 


n 


, r ,■*) 


u„(r lj )0' 


28. 


where the O v [j were introduced in Equation 1 1 . The spin-orbit correlations are only 
summed in Equation 22 because of the extra computational expense of evaluating 
powers of L • S that would occur if it were inserted in the symmetrized product of 
Equation 23. 

The central f c (r) and noncentral u p (r ) pair correlation functions reflect the 
influence of the two-body potential at short distances, while satisfying asymptotic 
boundary conditions of cluster separability. Reasonable functions are generated by 
minimizing the two-body cluster energy of a somewhat modified NN interaction; 
this results in a set of eight coupled, Schrodinger-like, differential equations cor- 
responding to linear combinations of the first eight operators of Vjj, with a number 
of embedded variational parameters (10). The f c (r) is small at short distances, to 
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reduce the contribution of the repulsive core of , and peaks at an intermediate 
distance corresponding to the maximum attraction of Vjj , as illustrated in Figure 2 
for several nuclei. For the s-shell nuclei, f c (r) falls off at larger distances to keep 
the system confined. For example, in 4 He, [/ c (r)] 3 ~ exp(— Kr)/r at large r, as 
shown by the dashed line in Figure 2, where k corresponds to an a t + p sepa- 
ration energy. The noncentral u p (r) are all relatively small; the most important is 
the long-range tensor-isospin part u lT (r), also shown in Figure 2 for several nuclei, 
which is mainly induced by the one-pion-exchange part of i>, ; . 

The f'j k , f/j k , and t/,j* are three-nucleon correlations also induced by u, ; - (11). 
The U[ ]k ' are three-body correlations induced by the NNN interaction, which have 
the form suggested by perturbation theory: 

U i jk — e - x V‘J k (?ij ’fjk, ?ki ) , 29 . 

X 

with r = yr, y a scaling parameter, and e x a (small negative) strength parameter. 
A wave function somewhat simpler to evaluate than 'T v is given by 


I'J'r) = 


1 + E ”1 


TNI 

Uk 


i<j<k 


I^P>. 


30. 


where U j^ 1 is a slightly truncated TNI correlation. The Ty gives about 1-2 MeV 
less binding than the full 'Tv but costs less than half as much to construct, making 
it a more efficient starting point for the GFMC calculation (18) . 

The Jastrow wave function for the light p-shell nuclei is more complicated, 
as a number of nucleons must be placed in the unfilled p-shell. The LS coupling 
scheme is used to obtain the desired JM value of a given state, as suggested by the 
shell-model studies of light p-shell nuclei (5). Different possible LS combinations 
lead to multiple components in the Jastrow wave function. The possibility that the 
central correlations f c (rjj) could depend on the shells (s or p) occupied by the 
particles and on the LS coupling is also allowed for: 


l*7> =A 


r~I fijk 0 fss( r ij ) ]~ [ fsp ( r kl ) 

i<j<k i<j< 4 k<4<l<A 


E (&SM n /p/ [ " 1 ('-/m)|4>A(L5[n]yMrr 3 ) 12 34:56...A) 
LSIn] ^ 4 <l<m<A 


31. 


The operator A indicates an antisymmetric sum over all possible partitions of the 
A particles into 4 s-shell and (A — 4) p-shell ones. The central correlation f ss (r ) 
is the f c (r) from the 4 He wave function. The f sp (r), shown in Figure 2 for Li 
nuclei, is similar to the f c (r) at short range, but with a long-range tail going to 
unity; this helps the wave function factorize to a cluster structure like a + d in 6 Li 
or a + 1 in 7 Li at large cluster separations. The f pp ^'\r) is similar to the deuteron 
(triton) fc(r) in the case of fi Li ( 7 Li). 
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The LS components of the single -particle wave function are given by 

\<t> A (LS[n]JMTT 3 ) 12 34 : 56 ...A) = |<M0000) 123 4> n 

4 <1<A 


X 


n 

L 4 <1<A 


¥ 1 m/i^al) 


X 

>< 

E 

i 

[\ m °) 

) 

- LM L [n] 

L 4</<A 

V z / J 

SMs J 


JM 



32. 


The (pp^iRak) are p-wave solutions of a particle in an effective a + N potential 
and are functions of the distance between the center of mass of the alpha core and 
nucleon k\ they may be different for different LS[n] components. A Woods-Saxon 
potential well. 


V aN (r) = V p LS 


1 + exp 




33. 


where V p s , R p , and a p are variational parameters, and a Coulomb term if appro- 
priate are used. 

For A > 7 nuclei, with three or more p-shell particles there are multiple ways 
of coupling the orbital angular momentum, L, spin, S, and isospin, T, to obtain 
the total (J n \T) of the system. The spatial permutation symmetry, denoted by the 
Young pattern [«], is used to enumerate the different possible terms, as discussed 
in Appendix 1C of Bohr & Mottelson (37). The different possible contributions to 
A = 6-8 nuclei are given in Table 1, with the corresponding spin states of highest 
spatial symmetry for each nucleus. 


TABLE 1 Permutation symmetry terms for LS-coupling in A = 6-8 nuclei and 
corresponding spin states 


A 

[//] 

L 

(T,S) 

Highest symmetry states 

6 

[2] 

0,2 

(L 0)(0, 1) 

6 He(0+ 2+), 6 Li(l+ 2+, 3+) 


[11] 

1 

(1. D(0. 0) 

6 He(l+) 

7 

[3] 

1,3 

(id) 

’Li(L. f-. I-) 


[21] 

1, 2 

( 2 ’ 2 )( 2 ’ 2 )( 2 ’ 2 ) 

!-. i“) 


[111] 

0 

(i i)(i> 0 


8 

[4] 

0,2,4 

(0, 0) 

8 Be(0+ 2+, 4+) 


[31] 

1,2,3 

(1, 1)(L0)(0, 1) 

8 Li(0 + - 4 + ), 8 Be(l+, 3+) 


[22] 

0,2 

(2, 0)(1, 1)(0, 2)(0, 0) 

8 He(0+, 2+) 


[211] 

1 

(2, 1)(L2)(1, 1 )( 1, 1)(1, 0)(0, 1) 

8 He(l+) 
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After other parameters in the trial function have been optimized, a series of 
calculations are made in which the j J >i.s\n\ of Equation 32 may be different in the left- 
and right-hand-side wave functions to obtain the diagonal and off-diagonal matrix 
elements of the Hamiltonian and the corresponding normalizations and overlaps. 
The resulting NxN matrices are diagonalized to find the Pls[»] eigenvectors, using 
generalized eigenvalue routines because the correlated 4V are not orthogonal. This 
allows us to project out not only the ground states but also excited states of the 
same (J n ;T) quantum numbers. For example, the ground states of 6 Li, 7 Li, and 
8 Li are obtained from 3 x 3, 5 x 5, and 7x7 diagonalizations, respectively. 

The energy expectation value of Equation 21 is evaluated using Monte Carlo 
integration. A detailed technical description of the methods used can be found 
elsewhere (10, 38). Monte Carlo sampling is done both in configuration space and 
in the order of operators in the symmetrized product of Equation 23 by following 
a Metropolis random walk. The expectation value for an operator O is computed 
with the expression 

= E P .JdR[^l(R)0%(R)/W pg (R)]W pg (R) ^ 

where samples are drawn from a probability distribution, W pq { R). The subscripts p 
and q specify the order of operators in the left- and right-hand-side wave functions, 
and the integration runs over the particle coordinates R = (rj, r 2 , . . . , r^). The 
probability distribution is constructed from the T P of Equation 23: 

W pq (R) = |Re«*; p (R)'k P , 3 (R)))|. 35. 

This is much less expensive to compute than using the full wave function of 
Equation 22 with its spin-orbit and operator-dependent three -body correlations, 
but it typically has a norm within l%-2% of the full wave function. 

Expectation values have a statistical error that can be estimated by the standard 
deviation a : 


cr 


\{ 0 2 )-{ 0 ) 


21 1/2 


N - 1 


36. 


where N is the number of statistically independent samples. Block averaging 
schemes can also be used to estimate the autocorrelation times and determine 
the statistical error. 

The wave function T can be represented by an array of 2 A x ( 7 ) complex 
numbers. 


$(R) = ^fc(R)|«), 37. 

a 

where the i//«(R) are the complex coefficients of each state |a) with specific third 
components of spin and isospin. This gives vectors with 96, 1280, 4480, and 14,336 
complex numbers for 4 He, 6 Li, 7 Li, and 8 Li, respectively. The spin, isospin, and 
tensor operators contained in the two-body correlation operator U , s and in 
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the Hamiltonian are sparse matrices in this basis. For forces that are largely charge- 
independent, as is the case here, this charge-conserving basis can be replaced with 
an isospin-conserving basis that has N(A, T) — 2 A x /(A, T) components, where 


HA, T) = 


2T+1 / A \ 

\a + t + A\a + t)- 


38. 


This reduces the number of vector elements to 32, 320, 1792, and 7168 for the 
cases given above — a significant savings. In practice, isospin operators are more 
expensive to evaluate in this basis, but the overall savings in computation is still 
large. In addition, the even-A nuclei can be computed in the M — 0 state, which 
requires only half the spin components; the others can be found by time-reversal 
invariance. 

Expectation values of the kinetic energy and spin-orbit potential require the 
computation of first derivatives and diagonal second derivatives of the wave func- 
tion. These are obtained by evaluating the wave function at 6A slightly shifted 
positions of the coordinates R and taking finite differences (as discussed in 10). 
Potential terms quadratic in L require mixed second derivatives, which can be 
obtained by additional wave-function evaluations and finite differences. 

The rapid growth in the size of ^ as A increases is the chief limitation, both 
in computer time and memory requirements, on extending the current method 
to larger nuclei. The cost of an energy calculation for a given configuration also 
increases as the number of pairs, P = A( A — 1), because of the number of 
pair operations required to construct T, and roughly as the number of particles, 
A, because of the number of wave -function evaluations required to compute the 
kinetic energy by finite differences. Table 2 shows these factors for a number of 
nuclei, with a final overall product showing the cost relative to that for a 8 Be 
calculation. Some initial calculations for A = 9, 10 have been made, but these are 
at the limit of current computer resources, so although 12 C should be reached in 
a few years, it may be the practical limit for this approach. Calculations of nuclei 
such as ie O or 40 Ca will require other methods, such as auxiliary-field diffusion 
Monte Carlo (AFDMC) which samples the spin and isospin components of the 
wave function by introducing auxiliary fields, as the plain Monte Carlo samples the 
spatial part. Also shown in Table 2 is the cost for several pure neutron problems, 
such as an eight-body drop (with external confining potential) or 14 or 38 neutrons 
in a box (with periodic boundary conditions) as neutron matter simulations. The 8 n 
drop has been calculated (39, 40) by VMC and GFMC, and initial box simulations 
have been done for 14 neutrons by GFMC (J. Carlson, private communication) 
and for up to 54 neutrons by AFDMC (20). 

A major problem arises in minimizing the variational energy for p-shell nuclei 
using the above wave functions: There is no variational minimum that gives reason- 
able rms radii. For example, the variational energy for 6 Li is slightly more bound 
than for 4 He, but is not more bound than for separated 4 He and 2 H nuclei, so the 
wave function is not stable against breakup into a + d subclusters. Consequently, 
the energy can be lowered toward the sum of 4 He and 2 H energies by making the 
wave function more and more diffuse. Such a diffuse wave function would not be 
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TABLE 2 Approximate scaling of VMC and 
GFMC calculations with system size* 



A 

P 

2 A x / (A,T) 

n/nc 8 Be) 

3 h 

3 

3 

8x2 

0.0004 

4 He 

4 

6 

16 x 2 

0.001 

5 He 

5 

10 

32 x 5 

0.020 

6 Li 

6 

15 

64 x 5 

0.036 

7 Li 

7 

21 

128 x 14 

0.66 

8 Be 

8 

28 

256 x 14 

1 

9 Be 

9 

36 

512x42 

17 

10 B 

10 

45 

1024 x 42 

24 

12 C 

12 

66 

4096 x 132 

530 

16q 

16 

120 

65536 x 1430 

2 x 10 5 

40 Ca 

40 

780 

7 x 10 21 

3 x 10 2 ° 

8 n 

8 

28 

256 x 1 

0.071 

14 n 

14 

91 

16384 x 1 

26 

38 n 

38 

703 

3 x 10“ 

9 x 10 9 


*The number of pairs is P; the state-vector size is the product of 
the number of spin states and number of isospin states, I(A,T). As 
discussed in the text, this number may be halved for even-A nuclei. 
The last column is the relative difficulty scaled to 8 Be, assuming 
M = 0 states are computed for even-A nuclei. 


useful for computing other nuclear properties, or as a starting point for the GFMC 
calculation, so the search for variational parameters is constrained by requiring the 
resulting point proton rms radius, r p , to be close to the experimental values for 6 Li 
and 7 Li ground states. For other A = 6-8 ground states, and all the excited states, 
the trial functions contain minimal changes to the 6 Li and 7 Li wave functions, with 
the added requirement that excited states should not have smaller radii than the 
ground states. The final step is always the diagonalization of the Flamiltonian in 
the Pls[h] mixing parameters. 

4. GREEN’S FUNCTION MONTE CARLO 

The GFMC method (3,41) projects out the exact lowest-energy state, To, for a 
given set of quantum numbers, using To = lim^oc exp [— (H — E’ojr] Ty, where 
Ty is an initial trial function. If the maximum r actually used is large enough, the 
eigenvalue Eq is calculated exactly, whereas other expectation values are generally 
calculated neglecting terms of order |To — Tyj 2 and higher (18). In contrast, the 
error in the variational energy, Ey, is of order |To — Tyj 2 , and other expectation 
values calculated with Tj- have errors of order |Tq — Tjj. In the following, we 
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present a brief overview of modem nuclear GFMC methods; more detail may be 
found elsewhere (18, 19). 

We start with the Ty of Equation 30 and define the propagated wave function 

*(r), 

q/(r) = e - (H - Eo)T ^ T = q/j., 39. 

where we have introduced a small time step, r = n Ar; obviously 'Lit = 0) = Ty 
and Tit — > oo) = 'Lq. The Ttr) is represented by a vector function of R using 
Equation 37, and the Green’s function, G a p (R, R'), is a matrix function of R and 
R' in spin-isospin space, defined as 

G a/J (R,R') = (R,a\e- (H - Eo)AT \R',P). 40. 

It is calculated with leading errors of order (Ar) 3 . Omitting the spin-isospin indices 
a and ft for brevity, 4TR, , r) is given by 


v I / (R„,r) = y G(R„, R„_i) • • ■ G(Ri, Ro)'I / 7-(Ro)a , R„_i • • • c/Rit/Ro, 41. 

with the integral being evaluated stochastically. 

The short-time propagator should allow as large a time step Ar as possible, 
because the total computational time for propagation is proportional to 1/Ar. 
Earlier calculations (3, 17,41) used the propagator obtained from the Feynman 
formulae, in which the kinetic and potential energy terms of H are separately 
exponentiated. The main error in this approximation to G a ^ comes from terms in 
e~ HAT having multiple v,j, such as Vjj K u i; (Ar) 3 , where K is the kinetic energy, 
which can become large when particles i and j are very close owing to the large 
repulsive core in u, ; - . This requires a rather small Ar ~ 0.1 GeV” 1 . 

It has been found in studies of bulk helium atoms (42) that including the exact 
two-body propagator allows much larger time steps. This short-time propagator is 


G ff/j (R,R') = G 0 (R,R')M 




Sijtoj’ r h > 


, 8o,ij(r u , r' ( ) 


\P), 


42. 


where 

G 0 (R,R') = (R|e _A:Ar |R') 

gij is the exact two-body propagator, 

8ij(r u ,r' ij ) = {r u \e 


/ m 

3A 

pvn 

-(R-R) 2 ' 

_V 2jtTi 2 At _ 

CAL/ 

2Ti 2 At / m 


~ H " AT K), 


43. 


44. 


and gojj is the free two-body propagator. All terms containing any number of the 
same v,j and K are treated exactly in this propagator, as we have included the 
imaginary-time equivalent of the full two-body scattering amplitude. It still has er- 
rors of order (Ar) 3 , but they are from commutators of terms such as V/j K %(Ar) 3 , 
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which become large only when both pairs ij and ik are close. Because this is a rare 
occurrence, a five-times-larger time step, Ar ~ 0.5 GeV - 1 , can be used (18). 

Including the three -body forces and the Eq in Equation 40, the complete prop- 
agator is given by 


G af i(R, R') = e EoAT G Q (R, R')exp 


£( V &( R > + V &( R ')) 


Ar 

~Y 


(a|/ 3 (R)|y>(y| 




KJ 


/ 8o,ij(T tj , r' ) 


|5><5|/ 3 (R')I>S), 


45. 


where 


h(R) = 


At x -, 


46. 


and V? jk V% + V% represents, in general, all noncentral VJ/* terms. With the 
exponential of V/J k expanded to first order in Ar, there are additional error terms 
of the form V^ k V^y k , (Ar) 2 . However, they have negligible effect because V* k has 
a magnitude of only a few MeV. 

Quantities of interest are evaluated in terms of a “mixed” expectation value 
between 4/ and 'l'(r): 

(W r | 0 |W(r)> 

1 )mixed_ <4 r |4(r)) 


f dP „ 4.}. (R„ ) O G (R„ , R„_i) • • ■ G(Ri, Ro)'lQ(Ro) 
/ r/P„'I , y(R„)G(R„, R„_i) • • • G(Ri, Rol'lGlRo) 


where P„ = Ro, Ri, . . . , R„ denotes the path, and dV n = c/Rot/Ri , . . . , d R„, with 
the integral over the paths being carried out stochastically. The desired expectation 
values would, of course, have 4(r) on both sides; by writing 4(r) = 4 r + 
A vp ( r ) and neglecting terms of order [A (r j] 2 , we obtain the approximate expre- 
ssion 


<0(r)> = 


MW) 

<4(t)|4(t)> 


(G(r))mixed + [(G( T )>mixed — (O)t], 


48. 


where (O)t is the variational expectation value. More accurate evaluations of 
(G(r)} are possible (43), essentially by measuring the observable at the mid-point 
of the path. However, such estimates require a propagation twice as long as the 
mixed estimate and require separate propagations for every (O) to be evaluated. 
The nuclear calculations published to date use the approximation of Equation 48. 

A special case is the expectation value of the Hamiltonian. The (H(r)) m ; xe d can 
be re-expressed as (44) 


(^f(Q)mixed 


(\\J T \e~ (H ~ E °' ,T ^ 2 He ~ (H ~ Eo)T ^ 2 1 4 7 - } 

(TQ \ e -lH-E 0 )T/2 e ^H-E 0 )T/2 | vpj.) - E °’ 


49. 
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since the propagator exp[— (H — Eq)t] commutes with the Hamiltonian. Thus 
( // ( t) ) m ixed approaches Eq in the limit r -> oo, and the expectation value obeys 
the variational principle for all r. 

The AVI 8 interaction contains terms that are quadratic in orbital angular mo- 
mentum, L. These terms are, in essence, state- and position-dependent modifica- 
tions of the mass of the nucleons. If they are included in the calculation of g tJ , 
then the ratio gij/go,ij in Equation 43 becomes unbounded for large r, ; - or rL 
and the Monte Carlo statistical error will also be unbounded. Hence the GFMC 
propagator is constructed with a simpler isoscalar interaction, H' , with a i , ' / that 
has only eight operator terms, [1, er, cr S,j, L • S]®[1, r,- -r ; ], chosen such that 
it equals the Cl part of the full interaction in all S- and P-waves and in the deuteron. 
The v'jj is a little more attractive than Vij, so a V- ]k adjusted so that (//'} ~ ( H ) 
is also used. This ensures the GFMC algorithm will not propagate to excessively 
large densities because of overbinding. Consequently, the upper bound property 
applies to (//'( r)), and ( H — H ') must be evaluated perturbatively. 

Another complication that arises in the GFMC algorithm is the “fermion sign 
problem.” This arises from the stochastic evaluation of the matrix elements in 
Equation 47. The G(R;, R,-_i) is a local operator and can mix in the boson so- 
lution. This has a (much) lower energy than the fermion solution and thus is 
exponentially amplified in subsequent propagations. In the final integration with 
the antisymmetric Ty, the desired fermionic part is projected out, but in the pres- 
ence of large statistical errors that grow exponentially with r. Because the number 
of pairs that can be exchanged grows with A, the sign problem also grows expo- 
nentially with increasing A. For A > 8, the errors grow so fast that convergence 
in r cannot be achieved. 

For simple scalar wave functions, the fermion sign problem can be controlled 
by not allowing the propagation to move across a node of the wave function. Such 
“fixed-node” GFMC provides an approximate solution that is the best possible 
variational wave function with the same nodal structure as Ty-. However, a more 
complicated solution is necessary for the spin- and isospin-dependent wave func- 
tions of nuclei. In the past few years a suitable “constrained path” approximation 
has been developed and extensively tested, first for condensed matter systems (45) 
and more recently for nuclei (19). The basic idea of the constrained-path method 
is to discard those configurations that, in future generations, will contribute only 
noise to expectation values. If the exact ground state To) were known, any con- 
figuration for which 

T(P„) + T 0 (R„) = 0, 50. 

where a sum over spin-isospin states is implied, could be discarded. Here P„ 
designates the complete path (Ro, . . . , R„) that has led to R„. The sum of these 
discarded configurations can be written as a state |T,y}, which obviously has zero 
overlap with the ground state. The T/ contains only excited states and should 
decay away as r -> oo; thus, discarding it is justified. Of course, the exact To is 
not known, and so configurations are discarded with a probability such that the 
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average overlap with the trial wave function vanishes, 

(’J'dl'Ey) = 0. 51. 

This procedure has been tested many times (19), and its results are usually con- 
sistent with unconstrained propagation, within statistical errors. However, a few 
cases have been found in which the constrained propagation converges to the 
wrong energy (either above or below the correct energy). Therefore, a small num- 
ber, n u — 10 to 20, of unconstrained steps are made before evaluating expectation 
values. These few unconstrained steps, out of typically 400 total steps, appear to be 
enough to damp out errors introduced by the constraint but do not greatly increase 
the statistical error. 

Figure 3 shows the progress with increasing r of typical constrained GFMC 
calculations, in this case for various states of s Be. The values shown at r = 0 are 
the VMC values using the 'Ey of Equation 30; VMC values with the better >E v 
of Equation 22 are about 1 .5 Me V below these. The GFMC very rapidly makes a 
large improvement on these energies; by r = 0.01 MeV” 1 , the 'Ey energies have 
been reduced by 6.5 to 8 MeV. This rapid improvement corresponds to the removal 
of small admixtures of states with excitation energies ~1 GeV from 'Ey. Averages 
over typically the last nine r values are used as the GFMC energy. The standard 
deviation, computed using block averaging, of all of the individual energies for 
these r values is used to compute the corresponding statistical error. The solid lines 
show these averages; the corresponding dashed lines show the statistical errors. 
The 0+ ground state, 1 + , and 3 + states of s Be have widths of less than 250 keV, 
and their E(r) appear to be converged and constant over the averaging regime. 
The 2 + state has a width of 1.5 MeV and the 4 + state’s width is 3.5 MeV; the 
E( r) for the 2 + state might be converged, but that of the 4 + state is clearly steadily 
decreasing. Reliable estimates of the energies of broad resonances require the use 
of scattering-state boundary conditions (see Section 7), which have not yet been 
implemented for more than five nucleons. 

Tests with different 'Ey as starting points for the GFMC calculation are described 
by Pudliner et al. (18). The most crucial aspect in choosing 'Ey is that it have the 
correct mix of spatial symmetries, i.e., the j J >i_s\n\ admixtures. This is because the 
limited propagation time is not sufficient for the GFMC algorithm to filter out 
low-lying excitations with the same quantum numbers. This issue is unimportant 
in 4 He, where the first 0 + excited state is near 20 MeV, but in 6 Li, the first 1 + excited 
state is at only 5.65 MeV; other p-shell nuclei have similar low-lying excitations. 
Otherwise, the GFMC algorithm is able to rapidly correct for some very poor 'Ey 
whose variational energies are actually positive. However, a good 'Ey helps to 
keep the error bars small at larger r . The most efficient balance between speed of 
construction of 'Ey and smallest number of samples needed to achieve a given error 
bar seems to be given by Equation 30. A number of other tests of the GFMC algo- 
rithm and its ability to determine nuclear radii are described elsewhere (18, 19). 

The calculations described here are computationally intensive and would not 
have been possible without the advent of parallel supercomputers. The initial 
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studies (17) of 6 Li required ~2000 node hours on an IBM SP1 to propagate 10,000 
configurations to r = 0.06 MeV -1 . Today, with improvements in the algorithm 
and after much effort to optimize the computer codes, the same calculation requires 
about 15 node hours on a third generation IBM SP. However, a present calculation 
for 8 Li, requiring 10,000 configurations to get a reasonable error bar, takes 1000 
node hours, so forefront computer resources remain essential for this program. 


5. ENERGY RESULTS 

There are a number of accurate many-body methods for evaluating the binding 
energy of the s-shell nuclei, 3 H, 3 He, and 4 He, using realistic nuclear forces. 
These include Faddeev in configuration space (46) (FadR), Faddeev and Faddeev- 
Yakubovsky in momentum space (47) (FadQ), and hyperspherical harmonics (48) 
(HH), correlated hyperspherical harmonics (49) (CHH) and pair-correlated hyper- 
spherical harmonics (50) (PHH). Table 3 shows results of these methods for the 
AVI 8 and AV18/UIX Hamiltonians. We observe that the VMC upper bounds are 
generally 1.5%-2.0% above the GFMC energies, which in turn are 0.25%-1.0% 
above the Faddeev and HH results. The disagreement between GFMC and the 
Faddeev and PHH methods may be attributable to the fact that GFMC is propa- 
gated with an H' and a small piece of the Hamiltonian is computed in perturbation. 
At this fine level of comparison, one must worry about how features such as T = | 
admixtures and n — p mass differences in the trinucleon ground state are treated 
in each calculation. 

Table 3 also shows the necessity of including (at least) a three-nucleon potential 
in the Hamiltonian in order to reproduce the A = 3,4 experimental binding energies. 
This is true for all the other modern NN potentials (46, 47). The local potentials 
AVI 8, Reid93, and Nijm II give very similar results, while the slightly nonlocal 
Nijm I gives l%-3% more binding, and the more nonlocal CD Bonn gives 5%-8% 
more binding, or about halfway between the local potential values and experiment. 
At present, no two- plus three-nucleon potential combination gives an exact fit to 
both trinucleon and 4 He energies, but several combinations, like AV18/UIX, come 


TABLE 3 Energies of 3 H and 4 He in MeV, computed by VMC and GFMC compared to 
momentum-space Faddeev (47), configuration-space Faddeev (46) and various hyperspherical 
harmonics methods (A. Kievsky, private communication) 


Hamiltonian 

Nucleus 

VMC 

GFMC 

FadQ 

FadR 

HH* 

Expt. 

AV18 

3 H 

—7.50(1) 

-7.61(1) 

-7.623 

-7.62 

-7.618 

8.482 


4 He 

-23.72(3) 

-24.07(4) 

-24.28 


-24.18 

28.296 

AV18/UIX 

3 H 

-8.32(1) 

-8.46(1) 

-8.478 


-8.475 

8.482 


4 He 

-27.78(3) 

-28.33(2) 

-28.50 


-28.1 

28.296 


*The HH column values are PHH for 3 H; HH for 4 He, AV18; and CHH for 4 He. AV18/UIX. 
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quite close. It also appears that it may be easier to fit 3 He and 4 He simultaneously, 
rather than 3 H and 4 He, perhaps because of the better data constraints on the pp 
interaction compared with the nn interaction (47). In principle, there should also be 
four-body forces, but their contribution must be quite small, and it is impossible 
to unambiguously identify a need for such terms until a more thorough survey 
of possible three-nucleon potentials is made and the discrepancies between the 
various many-body methods are resolved. 

Figure 4 compares VMC (using the full 4V of Equation 22) and GFMC calcu- 
lations for the AV18/UIX Hamiltonian and also shows experimental energies. All 
GFMC energies in this review are from Pieper et al. (22); VMC energies are from 
Pudliner et al. (18) and Wiringa et al. (19). All particle-stable or narrow-width 
(T < 150 keV) states for 4 < A < 8 are shown, except isobaric analogs. We 
see that the VMC calculations get progressively worse in the p-shell compared 
with the s-shell, being on the order of 10%-15% above the final GFMC results, 
for nuclei with N ~ Z. In absolute terms, 4V misses roughly an extra 1.5 MeV 
of binding for each p-shell nucleon that is added. The VMC results also fail to 
reproduce important qualitative features of the GFMC calculations; for example, 
6 Li is stable against breakup into a + d for this Hamiltonian, but the VMC calcu- 
lation shows it unbound by ~2 MeV (the black dashed lines show the indicated 
thresholds computed using the subcluster energies appropriate to each calculation 
or experiment). Clearly, there is some significant feature of light p-shell nuclei that 
is not yet included in the trial wave functions. 

Table 4 shows the GFMC energy values for several Hamiltonians. The table and 
Figure 4 show that the AV18/UIX Hamiltonian, which gives excellent energies, 


TABLE 4 GFMC and experimental energies of A = 3-8 ground 
states in MeV. Energies of artificially confined 7- and 8-neutron 
drops are also shown 



AV18 

AV18/UIX 

AV18/IL2 

Expt.* 

3 H(§+) 

-7.61(1) 

-8.46(1) 

-8.43(1) 

-8.48 

4 He(0 + ) 

-24.07(4) 

-28.33(2) 

-28.37(3) 

-28.30 

6 He(0 + ) 

-23.9(1) 

-28.1(1) 

-29.4(1) 

-29.27 

6 Li(l+) 

-26.9(1) 

-31.1(1) 

-32.3(1) 

-31.99 

7 He(|-) 

-21.2(2) 

-25.8(2) 

-29.2(3) 

-28.82 

7 Li(|~) 

-31.6(1) 

-37.8(1) 

-39.6(2) 

-39.24 

s He(0 + ) 

-21.6(2) 

-27.2(2) 

-31.3(3) 

-31.41 

8 Li(2+) 

-31.8(3) 

-38.0(2) 

-42.2(2) 

-41.28 

s Be(0 + ) 

-45.6(3) 

-54.4(2) 

-56.6(4) 

-56.50 

7 n(^-) 

-33.47(5) 

-33.2(1) 

-35.8(2) 


8 n(0 + ) 

-39.21(8) 

-37.8(1) 

-41.1(3) 



*The experimental values are from References 51-53. 
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TABLE 5 GFMC and experimental (53, 55) splittings of states with the 
same L and S but different J, in MeV 




L 

S 

AV18 

AV18/UIX 

AV18/IL2 

Expt. 

5 h 

1 - 3 - 

2 2 

1 

1 

2 

0.6(1) 

1.1(2) 

1.3(2) 

1.20 

6 Li 

2+ -3+ 

2 

1 

0.8(1) 

0.9(1) 

2.2(2) 

2.12 

7 Li 

1 - 3 - 

2 2 

1 

1 

2 

0.5(2) 

0.3(2) 

0.6(3) 

0.47 

7 Li 

5 - 7 — 

2 2 

3 

1 

2 

0.7(2) 

0.8(2) 

2.2(3) 

2.05 

8 Li 

1+-2+ 

1 

1 

0.2(4) 

0.6(2) 

L7(4) 

0.98 

7 n 

3 - 1 - 

2 2 

1 

1 

2 

1.65(7) 

1.5(1) 

2.8(3) 



compared with experiment, for the s-shell nuclei, significantly underbinds the 
p-shell nuclei. The Li isotopes are stable against breakup into subclusters with 
AV18/UIX, but progressively more underbound as A increases. There is also an 
isospin problem, in that the He isotopes are even further off from the experimental 
values and in fact do not show stability against breakup into subclusters. 

The new Illinois NNN potentials were constructed to solve the p-shell binding 
problems. Figure 5 and Table 4 compare GFMC calculations for AV18 with no 
Vjjk and AV18/IL2 to the same experimental energies shown in Figure 4. The 
AV18/IL2 Hamiltonian does a good job of reproducing these energies; the rms 
deviation from experiment for these levels is only 360 keV, whereas it is 2.3 MeV 
for AV18/UIX. The AVI 8 values with no Vijk show the large contribution that 
NNN potentials make to these binding energies; for A — 8, the IL2 increases the 
binding energy by more than 10 MeV. 

Table 5 shows some splittings for states that in a simple shell-model picture have 
the same L and S but different total J. In all cases, the lowest state for each / is used. 
Some of the states, particularly the 5 He states, have large experimental widths; cal- 
culations such as those described in Section 7 would be more reliable. The splittings 
are computed for the AVI 8, AV18/UIX, and AV18/IL2 Hamiltonians and are com- 
pared to experimental values. The AV18 with no Vj J<; significantly underpredicts 
all the splittings except the j — | doublet in 7 Li; however, in this case, the rel- 
atively large statistical error for the small splitting makes any conclusion difficult 
(the computed splittings are the difference of two independent GFMC calculations 
whose statistical errors must be added in quadrature). The AV18/UIX Hamiltonian 
substantially improves the 5 He splitting, and a VMC study showed that an earlier 
Urbana V tJ k significantly increased the computed splitting in 15 N, resulting in much 
better agreement with experiment (54). However, the other splittings in Table 5 
are still underpredicted with AV18/UIX. AV18/IL2 results in good predictions of 
all the splittings except for the 1 + — 2 + doublet in 8 Li, which is overpredicted. 

Table 6 provides a detailed breakdown of the GFMC ground-state energies for 
AV 1 8/IL2 into kinetic, NN and NNN interactions. Because of the extrapolation of 
the mixed expectation values (Equation 48), these components are not as accurate 
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TABLE 6 Kinetic, K, and potential energy contributions in MeV. Shown are the total 
NN and NNN potential energies, i>, ; and V ijk \ the electromagnetic, vj - , and one-pion, vfj , 
parts of v if, and the two- and three-pion parts of V ljk 



K 

v a 

v ijk 

*>5 

< 

v$ 

v$ 

H 

51 

-59 

-1.5 

0.04 

-45 

-3.0 

0.18(1) 

He 

115(1) 

-136(1) 

-8.4(1) 

0.86 

-105 

-16.3(1) 

0.63(1) 

He 

147(2) 

-171(2) 

-11.5(3) 

0.87 

-127(1) 

-20.3(4) 

-0.91(6) 

Li 

160(2) 

— 187(2) 

-11.1(3) 

1.73 

-150(1) 

-19.8(4) 

-0.44(5) 

He 

175(3) 

— 199(3) 

-16.3(4) 

0.89 

-145(2) 

-25(1) 

-3.1(1) 

Li 

199(3) 

-232(3) 

-14.5(4) 

1.80 

-178(2) 

-25.6(6) 

-1.1(1) 

He 

190(3) 

—218(3) 

-16.3(5) 

0.89 

-153(1) 

-25.6(6) 

-4.0(1) 

Li 

242(2) 

-278(2) 

-20.6(4) 

1.93 

-211(1) 

-34.2(5) 

-3.8(1) 

Be 

256(4) 

—303(3) 

-21(1) 

3.32 

-234(2) 

-38.5(9) 

-0.9(2) 

n 

105(1) 

-59(1) 

-3.6(3) 

0.07 

-10 

-0.1(1) 

-5.4(3) 

n 

122(1) 

-73(1) 

-3.0(3) 

0.09 

-12 

0.3(1) 

-5.9(4) 


as the total energy and do not add up to the full amount; indeed, the sum of the 
individual kinetic and potential energies differs from ( H ) by the same amount that 
the GFMC has improved the 'Ey energy. Nevertheless, they give a good idea about 
the relative size of the different terms involved. There is a big cancelation between 
the kinetic and two-body potential terms. Consequently, although V,yr is less than 
8% of v,j in magnitude, its expectation value is up to 50% of the net binding. 
The net effect of the IL2 VJ/*> defined as the difference of the binding energies 
computed without and with Vijk, is at most 30%, as can be seen in Table 4. This 
difference between expectation value and net effect is due to the large change that 
Vijk induces in (K + Vjj). 

Among the subcomponents of i> i; , the one-pion exchange dominates, being 
70%-80% of the total i> /; for A > 3 nuclei. Similarly, the two-pion exchange 
is the dominant component of V ljk . The three-pion exchange term is small and 
repulsive in s-shell nuclei but attractive in the p-shell nuclei, which have T — | 
triples. It is this term that results in the large improvement of the Illinois models 
over the Urbana models for p-shell nuclei; however, its contributions are always 
less than 15 % of the two-pion Vijk- Finally, we note that the electromagnetic vf- is 
dominated by the Coulomb interaction between protons, but about 17% (8%) of 
its total contribution comes from the magnetic moment and other terms in He (Li) 
isotopes. 

Energy differences among isobaric analog states are probes of the charge- 
independence-breaking parts of the Hamiltonian. The understanding of these 
“Nolen-Schiffer energies” has been a theoretical problem for 30 years (56). The 
energies for a given isomultiplet of states can be expanded as 



Annu. Rev. Nucl. Part. Sci. 2001.51:53-90. Downloaded from arjoumals.annualreviews.org 
by 193.219.53.38 on 09/26/05. For personal use only. 


QUANTUM MONTE CARLO CALCULATIONS OF LIGHT NUCLEI 


75 


TABLE 7 GFMC isovector and isotensor energy coefficients, a n (A,T) 
(Equation 52), computed using AV18/IL2, in keV* 


a„(A,T) 

j^CSB 

v C i (PP) 

v vJ< 

y CSB y CD 

Total 

Expt. 

ai(3,p 

14 

649(1) 

29 

64(0) 

757(1) 

764 

fli(6,l) 

16 

1091(5) 

18 

47(1) 

1172(6) 

1173 

a 2 (6,l) 


166(1) 

19 

107(13) 

293(13) 

223 

ai(7,i) 

22 

1447(6) 

40 

79(2) 

1588(7) 

1644 

fli(7,|) 

18 

1337(6) 

12 

52(1) 

1420(8) 

1373 

02(7, |) 


137(1) 

7 

36(6) 

180(7) 

175 

fli(8, 1) 

23 

1686(5) 

24 

76(1) 

1810(6) 

1770 

a 2 (8, 1) 


141(1) 

4 

-3(8) 

143(8) 

128 

fli(8, 2) 

18 

1528(7) 

17 

59(1) 

1622(8) 

1659 

a 2 (8, 2) 


136(1) 

6 

38(5) 

180(5) 

153 


' the contributions shown are kinetic energy, A' <SB ; Coulomb potential, Vci(pp)', other 
electromagnetic, v y ' R and nuclear potential, u CSB -+- i> CD . The experimental value for a-JX. 1) 
is extracted using 16.80 MeV as the excitation of the pure T = 1 isobaric analog state in 8 Be (58). 


EaAT;) = a '^ A ’ T )Qn(T, T z ), 52. 

n<2T 

where Q q = 1, Q\ = 72, and Qi = \ (3T 2 — T 2 ) are isoscalar, isovector, and 
isotensor terms (57). The isovector and isotensor coefficients a n (A . T), and various 
contributions to them, are given in Table 7. These were calculated as expectation 
values in the T z — —T GFMC wave functions for AV18/IL2. The individual terms 
are (a) the effect of the neutron-proton mass difference on the kinetic energy, 
K csb ; (b) the proton-proton Coulomb potential, including the AV18 form factor, 
i,’ci (pp)', (c) all other electromagnetic terms, such as vacuum polarization and 
magnetic -moment terms, v Y ' R \ and (cl) the strong-interaction contributions — u CSB , 
which contributes to the isotensor coefficients, and n CD , which contributes to the 
isovector coefficients. 

The 3 H- 3 He mass difference is 757 keV with the AV18/IL2 Hamiltonian, in 
excellent agreement with the experimental value of 764 keV. The bulk of the 
difference is the pp Coulomb energy, but 108 keV comes from the other terms. 
Because of its very long range, the Vci(pp) dominates the fli(A, T ) of heavier 
nuclei even more, and small errors in the rms radius of the nucleus can result in 
changes in the Vqi( pp) contribution that can totally mask the effects of the other, 
much smaller, terms. 

Two (2 + ;0+l) states occur very close together in the spectrum of x Be at 16.6 and 
16.9 MeV excitation; these isospin-mixed states come from blending the (2 + ;l) 
isobaric analog of the x Li ground state with the second (2 + ;0) excited state. There 
are also fairly close (1 + ;0,1) and (3 + ;0,l) pairs at slightly higher energies in the 
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TABLE 8 VMC isospin-mixing matrix elements for 8 Be, in 
keV. computed using AV18/UIX 


r 

j^CSB 

r C i (PP) 

v 7 -* 

V CSB 

Eoi 

Expt.* 

2+ 

2 

62 

19 

26 

109(4) 

149 

1 + 

1 

39 

0 

15 

55(2) 

120 

3+ 

1 

33 

15 

13 

62(2) 

63 


* Experimental values are from Reference 58. 


8 Be spectrum. The isospin-mixing matrix elements that connect these pairs of 
states, 


E 0 l (J)= ('T(/ + ;0)|7/|'T(/+;1)>, 53. 

have been computed using VMC wave functions for the AV18/UIX Hamiltonian 
(19). Table 8 shows results for the Eqi(J). The experimental values are determined 
from the observed decay widths and energies (58). The dominant contribution, 
from the Coulomb potential, typically accounts for less than half of the matrix 
element. We see that the remaining part of the electromagnetic interaction and the 
strong CSB interaction can provide a significant boost, although the experimental 
mixing is still underpredicted by ~20%. It appears that these mixing elements are 
a more sensitive test of the small CSB components of the Hamiltonian than are the 
isobaric analog energy differences. 

6. DENSITY AND MOMENTUM DISTRIBUTIONS 

The one- and two-nucleon distributions of light p-shell nuclei are interesting in a 
variety of experimental settings. For example, the Borromean 6 He and 8 He nuclei 
are popular candidates for study as “halo” nuclei, whose last neutrons are weakly 
bound. In addition, the polarization densities of 6 Li and 7 Li are important because 
of possible applications in polarized targets. 

Table 9 shows point proton rms radii and quadrupole moments, as computed by 
GFMC for the AV18/IL2 model, along with experimental values. The experimental 
charge radii have been conveited to point proton radii by removing the proton 
and neutron (r 2 ) of 0.743 and —0.116 fm 2 , respectively. The radii are in good 
agreement with experiment, at least partly because the AV 1 8/IL2 model reproduces 
the binding energies of these nuclei very well. The quadrupole moments have been 
calculated in impulse approximation; two-body charge contributions are expected 
to provide only a few percent correction. The agreement with experiment is again 
fairly good, with the exception of the 6 Li quadrupole moment, which involves 
a delicate cancelation between the contributions from the deuteron quadrupole 
moment and the D-wave part of the a + d relative wave function. In general, 
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TABLE 9 GFMC values for point proton rms radii 
(in fm) and quadrupole moments (in fm 2 ) in impulse 
approximation, computed with AV18/IL2 



<g 2 > 1/2 


Q 



GFMC 

Expt.* 

GFMC 

Expt.* 

3 h 

1.59(0) 

1.60 



3 He 

1.76(0) 

1.77 



4 He 

1.45(0) 

1.47 



6 He 

1.91(1) 




6 Li 

2.39(1) 

2.43 

-0.32(6) 

-0.083 

7 Li 

2.25(1) 

2.27 

-3.6(1) 

-4.06 

7 Be 

2.44(1) 


-6.1(1) 


8 He 

1.88(1) 




8 Li 

2.09(1) 


3.2(1) 

3.11(5) 

8 B 

2.45(1) 


6.4(1) 

6.8(2) 


* Experimental values are from References 55, 59-61. 


quadrupole moments are difficult to calculate accurately with quantum Monte 
Carlo methods because they are dominated by the long-range parts of the wave 
functions, which contribute very little to the total energy that VMC and GFMC 
are both designed to optimize. 

Table 10 shows calculated and experimental isoscalar and isovector magnetic 
moments, as defined with the convention used in Equation 52; thus, the isovector 
values for the T — \ cases are —2 times those often quoted in the literature. 
(For the A — 8, T— 1 nuclei, the experimental isoscalar and isovector moments 
are obtained from the sum and difference of the values for B and Li, since the 
magnetic moment of the T = 1 , J K = 2 + state in Be is not measured.) The 
moments were calculated as expectation values in the GFMC wave functions with 
T z — —T in impulse approximation. These impulse isoscalar moments are quite 
close to experiment, and it is expected that two-body current contributions are 
small (62). For isovector moments, however, there can be significant corrections 
from meson-exchange contributions. We note that the corrections needed for the 
A — 7,8 isovector moments are the same sign and magnitude as for A = 3, so an 
exchange-current model that fixes the s-shell (63) is likely to work for the light 
p-shell nuclei also. 

In Figure 6, we present proton and matter (proton+neutron) densities for the 
stable helium isotopes, as calculated with GFMC for AV18/IL2. Nucleon densities 
are calculated as simple 5-function expectation values, with possible spin and/or 
isospin projectors. For example, the proton density is given by 
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TABLE 10 GFMC values for isoscalar and isovector 
magnetic moments (in p N ) in impulse approximation using 
AV18/IL2 



T Isoscalar 


T Isovector 



GFMC 

Expt.* 

GFMC 

Expt.* 

3 He- 3 H 

1 

2 

0.403(0) 

0.426 

-4.330(1) 

-5.107 

6 Li 

0 

0.817(1) 

0.822 



7 Be- 7 Li 

1 

2 

0.894(1) 

0.929 

-3.93(1) 

-4.654 

8 B- 8 Li 

1 

1.276(1) 

1.345 

0.369(9) 

-0.309 


* Experimental values are from References 51-53, 61. 


Pp(r) = E i ~^p Ls(r - |r; - R ™DI' I '>- 54. 

i 

The blue symbols and curves in Figure 6 show results for 6 He and the red ones 
for x He; open symbols give the proton distributions, which are also interpreted 
as the alpha “core” density; full symbols are the total matter densities. It can be 
seen that, as more neutrons are added, the tails of the matter distributions broaden 
considerably because of the relatively weak binding of the p-shell neutrons. In 
addition, the central neutron and proton densities decrease rather dramatically. 
This effect does not necessarily require any changes to the alpha core but can be 
understood at least partially from the fact that the alpha no longer sits at the center 
of mass of the entire system. The motion relative to the center of mass spreads out 
the mass distribution relative to that of 4 He. 

Figure 6 also shows two attempts to extract the densities of 6x He from scattering 
of beams of these short-lived nuclei from proton targets (64). The analysis of Alk- 
hazov et al. (64) (dashed curves) used Glauber theory and assumed proton (core) 
and matter density distributions that were varied to fit the measured cross sections. 
The later work of Al-Khalili & Tostevin (65) (dot-dashed curves) improved on this 
by doing the Glauber theory using model wave functions that contained correla- 
tions between the valence neutrons. The GFMC calculations definitely prefer the 
latter analysis. 

The proton-proton distribution function is defined by 

Ppptr) = 4-2 W\ - l r '- 55 - 

i<j 

These distributions are directly related to the Coulomb sum measured in inclusive 
longitudinal electron scattering; such measurements in 3 He have been used to put 
constraints on the p pp (rjj), and realistic calculations agree with the experimental 
results (66). The behavior of p pp (r) at short distances is largely determined by the 
repulsive core of the NN potential and is nearly independent of the nucleus, but at 
larger distances it is determined by the size of the nucleus. 
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It is interesting to compare p pp for 4 He to that of 6 He and x Hc. Figure 7 shows 
these distribution functions, again calculated with GFMC for the AV 1 8/IL2 model. 
These nuclei each have just one pp pair, which presumably is in the alpha core 
of 6,8 He. Unlike the one-body densities, these distributions are not sensitive to 
center-of-mass effects, and thus if the alpha core of 6,8 He is not distorted by the 
surrounding neutrons, all three p pp distributions in the figure should be the same. 
We see that the pp distribution spreads out slightly with neutron number in the 
helium isotopes, with an increase of the pair rms radius of ~4% in going from 4 He 
to 6 He and of ~8% to s He. Although this could be interpreted as a swelling of the 
alpha core, it might also be caused by the charge-exchange (r,- ■ r f) correlations, 
which can transfer charge from the core to the valence nucleons. Since these 
correlations are rather long-ranged, they can have a significant effect on the pp 
distribution. VMC calculations of 4 He with wave functions modified to give p pp 
distributions close to those of 6,8 He suggest that the alpha cores of 6,8 He are excited 
by ~80 and ~350 keV, respectively. 

In general, VMC calculations give one -body densities very similar to the GFMC 
results, although the two-body densities may differ by up to 10% in the peak, 
with the GFMC having a somewhat sharper structure. Because the VMC wave 
functions are simpler and less expensive to construct, they have been used in a 
number of applications where the single-particle structure is dominant. All these 
calculations were carried out with wave functions for the AV18/UIX Hamiltonian, 
but we expect that GFMC calculations with the improved AV18/IL2 model will 
not qualitatively alter the results. (Remember that the VMC wave functions are 
constrained to reproduce available experimental rms radii.) 

VMC calculations of elastic and inelastic electromagnetic form factors for 6 Li 
are shown in Figure 8. These have been made in impulse approximation and with 
meson-exchange contributions to the charge and density operators (12). The com- 
parison with data for the elastic longitudinal form factor, F L (q 2 ), is excellent, as is 
the E2 transition to the 3 + first excited state. The meson-exchange contribution cor- 
rections are small, but they stand out noticeably in the first minimum, where they 
significantly improve the fit to data. The elastic transverse form factor, Fjiq 1 ), is 
good up to the first zero but is a little too large in the region of the second maximum. 
The Ml transition to the 0 + isobaric analog state is also reproduced reasonably 
well. This kind of quantitative agreement with data has not been achieved in the past 
with either shell-model or a + d cluster wave functions. In fact, shell-model cal- 
culations normally require the introduction of effective charges, typically adding 
a charge of ~0.5e to both the proton and neutron, to obtain reasonable transition 
rates (74). No effective charges are used in the VMC calculations. 

Whereas electron scattering from nuclei is primarily sensitive to the proton 
distributions in nuclei, n~ scattering is most sensitive to the neutron distribu- 
tions. Pion inelastic scattering at medium energy (80 < L| ah < 300 MeV) is 
dominated by strong absorption owing to excitation of the A resonance, and is 
well described in distorted-wave impulse approximation (DWIA). An analysis of 
meson-factory data for p-shell nuclei based on the Cohen-Kurath shell model (5) 
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6 (deg) 

Figure 9 Differential cross sections for 7 Li(7r, it') scattering at E n = 164 Me V. 
Data are from Reference 75. 


showed that reasonable agreement with data could be achieved if the quadrupole 
excitation operator was enhanced by a factor of ~ 1.75-2.5 (76). Figure 9 shows 
a calculation of the differential cross sections for 7 Li(7r, it’) scattering to the first 
three excited states using VMC transition densities as input, with no enhance- 
ment factors (14). Solid lines show the full results, which are in good agreement 
with the data; dashed lines give the contribution from protons only. Because of 
the strong isospin dependence of the ttN scattering f-rnatrix, this reaction can be 
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a stringent test of the neutron densities predicted by the quantum Monte Carlo 
calculations. 

The VMC wave functions for the AV18/UIX model have also been used to cal- 
culate single -nucleon momentum distributions in many nuclei, as well as a variety 
of cluster-cluster overlap wave functions, such as (dp\t), ( dd\a }, and (at/| 6 Li> 
(77,78). Recently, the overlaps ( 6 He(7 7r ) + /;>(£ 7 )| 7 Li} for all possible p-shell 
states in r ’Hc were studied (13). The spectroscopic factors obtained from these 
overlaps are 0.41 to the ground state of 6 He and 0.19 to the 2 + first excited state. 
These factors are significantly smaller than the predictions of the Cohen-Kurath 
shell model (79), which gives values of 0.59 and 0.40, respectively. The Cohen- 
Kurath shell model requires that the possible 6 He + p states sum to unity within 
the p-shell, whereas in the VMC calculation, the correlations in the wave function 
push significant strength to higher momenta that cannot be represented as a 6 He 
state plus p-wave proton. The VMC overlaps were used as input to a Coulomb 
DWIA analysis of recent 7 Life, e'/?) 6 He data taken at NIKHEF (13) and were 
found to give an excellent fit to the data, as shown in Figure 10. 

The VMC wave functions are based on one-body parts that have a shell-model 
structure, namely four nucleons in an alpha core coupled to (A — 4) one-body 
l — l wave functions. However, the low-lying states of 8 Be exhibit a rotational 
spectrum and are believed to be well approximated as two alphas rotating around 



p m [MeV/c] 

Figure 10 Experimental (13) momentum distributions for the transitions to the 
ground state (circles) and first excited state (squares) in the reaction 7 Li ( e,e'p ) 6 He, 
compared to Coulomb DWIA calculations with the VMC wave functions (solid lines) 
for AV18/UIX. 
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their common center of mass. It is possible to recover this picture from the VMC 
wave functions by a modified Monte Carlo density calculation (19). 

The standard Monte Carlo method for computing one-body densities, p(r), is 
to make a random walk that samples | Vk(ri ,t 2 , . . r,t)| 2 and to bin r| , r 2 , . . . . r,i 
for each configuration in the walk. The density is then proportional to the number 
of samples in each bin. In the case of a / = 0 nucleus, this “laboratory” density 
will necessarily be spherically symmetric. 

The intrinsic density in body-fixed coordinates can be approximated by com- 
puting the moment-of-inertia matrix, AT of the A positions for each configuration. 
The eigenvalues and eigenvectors of A4 are found and a rotation to those prin- 
cipal axes is made. The resulting r'j, r' 2 , . . . , r' A is then binned. The eigenvector 
with the largest eigenvalue is chosen as the z' axis. This procedure will not pro- 
duce a spherically symmetric distribution, even if there is no underlying deformed 
structure, because almost every random configuration will have principal axes of 
different lengths, and the rotation will always orient the longest principal axis in 
the z ' direction. However, no artificial structure is introduced (19). 

When the above procedure is applied to the three lowest 8 Be states, a dramatic 
intrinsic structure is revealed, as shown in Figure 1 1 for the ground state. The figure 
shows contours of constant density plotted in cylindrical coordinates. The left side 
of the figure shows that the standard (lab-frame) density calculation is spherically 



Figure 11 Calculated density contours of s Be in the lab frame (left) and the 
intrinsic frame (right), labeled with densities in fm -3 . 
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symmetric for the J — 0 ground state. The right side of the figure shows that the 
intrinsic density has two peaks, with the neck between them having only one third 
the peak density; we regard these as two alphas. This assignment is strengthened by 
making the same construction for the J = 2 + and 4 + states; although the laboratory 
densities for these states (in M — J states) are quite different, the intrinsic densities 
are, within statistical errors, the same as the J — 0 intrinsic density (19). 

If the 0+, 2 + , and 4 + states are generated by rotations of a common deformed 
structure, then their electromagnetic moments and transition strengths should all 
be related to the intrinsic moments, which can be computed by integrating over the 
projected body-fixed densities. This is explored by Wiringa et al. (19), who find a 
generally consistent picture. 


7. LOW-ENERGY NUCLEON-NUCLEUS SCATTERING 


The calculations presented so far have treated resonant nuclear states as if they 
were particle-stable; that is, the VMC trial wave functions decay exponentially 
at large distances and the GFMC propagation does not impose a scattering-state 
boundary condition. As was shown in Figure 3, this GFMC propagation converges 
for unstable states that are narrow, but for wide states, such as 8 Be(4 + ), the E( r) 
keeps decreasing for increasing r. Such states should be computed using scattering- 
wave boundary conditions. 

The 5 He(i , | ) resonances have been studied in VMC (80) and GFMC (81) 
using scattering-wave boundary conditions. The VMC calculations use a trial 
function that is a straightforward generalization of the T v given in Section 3. 
In it, the 4Q contains one p-wave single-particle function that goes to zero at 
a specified (large) radius, R„. The variational energy is minimized subject to 
this boundary condition. The phase shift at the resulting scattering energy, E s = 
L[ 5 Hei R„)] — /ft 4 He), is then obtained from 


tan(<5/) = 


ji(kR„) 

ni(kR„)’ 


56. 


where k — (2 / j , E s ) 1 ^ 2 /ft, fu, is the a + n reduced mass, and y) and n; are the 
spherical Bessel functions (in this case l — 1). It is possible to generalize this 
formulation by specifying a logarithmic derivative at R„ , rather than having the 
wave function go to zero there. This method requires the computation of the energy 
difference, E s . If this is done with two different VMC calculations, one for 5 He with 
the specified boundary condition and one for 4 He, then the statistical errors of 
these two calculations must be added in quadrature. Instead, one can directly 
evaluate the energy difference in a VMC random walk that is controlled by the 
5 He wave function. In practice, this gives E s with a smaller statistical error than 
either of the separate statistical errors (80). The GFMC calculations start with 
the scattering-state trial function and use a modified propagator that preserves 
the boundary condition (81). Both methods result in a phase shift and energy 
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corresponding to a given boundary condition. By repeating the calculations for 
different boundary conditions, the phase shift as a function of energy can be mapped 
out. 

The VMC calculations of Carlson et al. (80) used older two- and three-nucleon 
potentials than have been used for the other calculations reported in this review. 
They reproduced the qualitative features of the experimental a + n \ and f 
phase shifts, but, in particular, gave a spin-orbit splitting about 1 MeV too small. 
The GFMC calculations (81) used a different but still somewhat old Hamiltonian 
and obtained values that are closer to experiment. These calculations need to be 
repeated with the modern Hamiltonians described in this review, and for other 
systems such aso' + rf, 6 He + n, and a + a. 

8. ASTRO PHYSICAL ELECTROWEAK REACTIONS 

Most of the key nuclear reactions in primordial nucleosynthesis (82) and solar neu- 
trino production (83) involve only the A < 8 nuclei. Many of these are electroweak 
capture reactions that are difficult or impossible to measure in the laboratory. With 
the high-precision few-body methods, such as PHH and Faddeev, or new effective 
field theory treatments, the reactions to s-shell final states have been calculated with 
unprecedented accuracy in the past few years. These include the 'H(p, e + v e ) 2 H 
(84) and 3 He(p, e + v e ) 4 He (85) weak capture reactions, as well as the 'H(n, y) 2 H 
(86), 2 H (n, yj’H, and 2 H ( /; . y) 3 He (87) radiative capture reactions. 

Recently, the VMC method has been applied to some of the radiative captures 
to p-shell final states, including 2 H(a, y) 6 Li (15), 3 H(a, y) 7 Li, and 3 He(a, y) 7 Be 
(16). As discussed in the previous section, many -nucleon scattering states can be 
constructed within the quantum Monte Carlo framework, but this has not yet been 
done for A > 6 nuclei. Thus, for these first studies, the appropriate electromagnetic 
matrix elements (primarily El and E2) are evaluated between a VMC six- or 
seven-nucleon final state and a correlated cluster-cluster scattering state that is not 
variationally improved. 

Because much of the capture takes place at long range, an important ingredient 
in these calculations is an asymptotically correct description of the six- or seven- 
body final state. This feature has to be implemented in the initial variational wave 
function because, as remarked above, the VMC and GFMC methods find their best 
wave functions by optimizing the energy, which is not sensitive to the long-range 
behavior. The asymptotic behavior can be imposed by requiring the one-body wave 
function in Equation 33 to satisfy the condition 

[cpp Sln] (r -> oo)]” oc W km (2yr)/r, 57. 

where Wi, :m (2yr ) is the Whittaker function for bound-state wave functions in a 
Coulomb potential and n is the number of p-shell nucleons. Here y 2 — 2/x 4 „ B± n /tr, 
with /j, 4„ and B 4 ,, the appropriate two-cluster reduced mass and binding energy. 
For 6 Li, B 42 — 1.47 MeV; for 7 Li or 7 Be, capture can be to the | ground or 4 
excited state with corresponding values for B43. 
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The initial-state wave functions are taken as elastic-scattering states of the form 


\f ar :LSJM) = A 


4>£(r aT )Y L M L (?aT) n G u\M"r ,! 


58. 


LSJM 


where curly braces indicate angular momentum coupling, A antisymmetrizes be- 
tween clusters, i j/ a is the 4 He ground state, and 1 //“' is the deuteron or trinucleon 
ground state in spin orientation ms- The G, 7 are identity operators if the nucleons 
i and j are in the same cluster; otherwise, they are a set of pair correlation oper- 
ators, including both central and spin-isospin-dependent terms, which introduce 
distortions in each cluster under the influence of individual nucleons from the other 
cluster. They are similar to the correlations discussed in Section 3, except that they 
revert to the identity operator at pair separations beyond about 2 fm. The corre- 
lations </> 77 are generated from optical potentials that describe the experimental 
phase shifts in cluster-cluster scattering (see 15, 16 for details). Because the G, ; 
go to unity, the \jj- aT has the same phase shifts as those generated by the optical 
potential. 

Evaluation of electromagnetic matrix elements between the initial scattering and 
final bound states can be split into two parts by noting that all of the energy depen- 
dence is contained in the relative wave function </r /4 and the transition operators. 
Thus, using a technique first applied in a VMC calculation of d + d radiative cap- 
ture (88), the matrix element for a given scattering partial wave can be written as 


„ LSJiJf 


,J/ (<7) = 


/»oo 

/ dx x 2 </> 7, r L (x) 

Jo 


x{\ls A fmf \T a (q)A 


S(x - r aT )Y™ L ( f aT )) n G^IVuC 


59. 


LSJ/Mi 


where 7); denotes any of the standard E f or M, operators. The integration over all 
coordinates except x can be calculated just once for each partial wave by Monte 
Carlo sampling, and the result can then be used to compute the full integral for as 
many energies as desired by recomputing (p^ ix) only. 

The result for the 2 H(qi, y ) 6 Li capture 5-factor is shown in Figure 12, where it 
is compared to a collection of direct and indirect data. The 5-factor falls rapidly 
at low energy because a pseudo-orthogonality between the initial and final states 
suppresses normal S-wave capture. Whenever there is such a suppression, it is im- 
portant to consider higher-order terms that might contribute, including relativistic 
corrections and two-body charge and current operators (4). This is possible with 
the VMC calculation because of the fully correlated ,4 -body wave functions that 
are used. In this case, a relativistic center-of-energy correction leads to an El con- 
tribution that dominates at low energy, whereas the resonance region and above is 
primarily E2. Results for ] H(u, y) 7 Li and 3 He(a, y ) 7 Be can be found elsewhere 
(16). 

This kind of VMC calculation is only a first investigation into the p-shell as- 
trophysical reactions. It should be improved in the future by using the more exact 
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GFMC wave functions for both bound and scattering states, and the Hamiltonians 
with improved NNN potentials. 


9. NEUTRON DROPS 

Neutron-rich nuclei are interesting both because of their importance in various 
astrophysical contexts, such as the r-process or neutron-star crusts, and because of 
the current interest in experiments with radioactive beams. These nuclei are often 
studied within the framework of mean- field models using Skyrme or other potential 
models. The parameters of these models are fit to experimentally known binding 
energies, that is, for situations with N ~ Z. In particular, the isospin dependence 
of the spin-orbit component of such potentials is considered to be not strongly 
constrained. Neutron drops offer the possibility of theoretical guidance for the 
isotopic dependence of such parameters. 

Neutron drops are systems of interacting neutrons confined in an artificial exter- 
nal well. Eight neutrons form a closed shell and single -particle spin-orbit splittings 
can be studied in drops of seven or nine neutrons. Pairing energies can be studied 
if six-neutron drops are also computed. Calculations of systems of seven and eight 
neutrons interacting with AV18/UIX were used as a basis for comparing Skyrme 
models of neutron-rich systems with microscopic calculations based on realistic 
interactions (39). The external one-body well used is a Woods-Saxon, 


Vi(r) = Y V ° ; 60. 

i 1 + exp[— (r,- - r 0 )/a 0 ] 

the parameters are Vo = —20 MeV, r 0 = 3.0 fm, and a () — 0.65 fm. Neither 
the external well nor the total internal potential (u, 7 + Vijk) is individually at- 
tractive enough to produce bound states of seven or eight neutrons; however, the 
combination does produce binding. 

Tables 4—6 show results for the neutron drops with this external well. The 
T = | nature of the S' T term of V 3n results in large contributions in the neutron 
drops. As a result, the seven-neutron drops computed with AV18/IL2 have double 
the spin-orbit splitting predicted by AV18/UIX. Thus, the spin-orbit splitting in 
neutron-rich systems depends strongly on the Hamiltonian used. 


10. CONCLUSIONS AND OUTLOOK 

Quantum Monte Carlo methods can now be used to obtain accurate (within 2%) 
energies for light p-shell nuclei up to A = 8. Initial calculations of A — 9, 10 
nuclei have been made by the authors, and studies of A = 11, 12 nuclei should 
be feasible in the next few years by variational and Green’s function Monte Carlo 
methods. This progress in the nuclear many -body problem is due both to the 
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rapid growth of computational power and the continuing evolution of algorithms. 
The new auxiliary-field diffusion Monte Carlo method, which samples spins and 
isospins by auxiliary fields and samples space by standard diffusion Monte Carlo, 
may be the key to computing the properties of even larger nuclear systems (20). 

Studies of the p-shell nuclei allow us to test nuclear forces in new ways not 
accessible in s-shell nuclei, particularly the odd partial waves of NN scattering 
and the T — | triples for NNN forces. They also give us many more cases in 
which to examine charge-dependent and charge-symmetry-breaking interactions. 
With these calculations of nuclear spectra, we see for the first time that nuclear 
structure, including both single-particle and clustering aspects, really can be ex- 
plained directly in terms of bare nuclear forces that fit NN data. A crucial ingre- 
dient for quantitative agreement is the addition of realistic NNN forces, including 
at least two terms beyond the standard long-range two-pion-exchange potential 
( 22 ). 

Many aspects of nuclear structure and reactions are described quantitatively in 
these studies. The energy differences between isobaric analog states are explained 
well once the complete electromagnetic and strong charge-independence-breaking 
forces deduced from NN scattering are included. Charge radii and quadrupole mo- 
ments agree with experiment, and we expect magnetic moments will also agree 
once the important two-body exchange currents are included, as they have been 
in s-shell nuclei (63). Elastic and transition form factors measured in electron 
scattering, and transition densities that are tested in pion scattering, are accu- 
rately predicted without the introduction of effective charges. Spectroscopic fac- 
tors are naturally quenched by the correlations in ,4 -body wave functions and 
we see that those correlations build up intrinsic cluster structure, as in the case 
of 8 Be. 

Initial calculations of nucleon-nucleus scattering states and electroweak capture 
reactions are encouraging. These two related problems will be a major activity for 
the quantum Monte Carlo studies over the next several years, with studies of 
scattering states in nuclei such as 7 He and 9 He, and of the cluster-cluster scattering 
states that enter into the astrophysical capture reactions. Additional reactions like 
7 Bet p, y f'ti and even 8 Be(a, y) 12 C should become feasible. Still other projects 
will include the study of weak decays and, going beyond the p-shell, the study of 
the unnatural-parity intruder states, which start to dip into the low-energy spectrum 
with the A = 9 nuclei and can be particle-stable by A — 10. 

Refining the nuclear Hamiltonian will remain a major aspect of future work. 
Extensive new NN scattering data taken since the Nijmegen partial-wave analysis 
have increased the data base to more than 6000 points, and a new CD Bonn 2000 
potential has been constructed that again achieves a / 2 /datum % j (92). Also, 
intensive searches are under way for NNN force signatures in N + d scattering 
experiments (93, 94), which may well indicate a need for additional terms beyond 
those contained in the Illinois models. Our ability to test new force models, as they 
become available, in light nuclei by accurate quantum Monte Carlo methods will 
continue to be an important tool for nuclear physics. 
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r(fm) 

Figure 2 Correlation functions in light nuclei: the central pair correlation, f c , and 
tensor-isospin correlation, u tT , for s-shell nuclei are shown with solid lines. The dashed 
lines show f;?( 3 H) and f 3 ( 4 He) to illustrate the long-range separation behavior. The dot- 
dashed lines show the central correlation between s- and p-shell nucleons, f sp , the u tr , 
and the one-body wave function, </>p SIn] , for 6 Li. 



Figure 3 GFMC E(r) using AV18/UIX for various states of s Be. Averages 
with their standard errors are shown by the solid and dashed lines. 
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Figure 4 VMC and GFMC energies using AV18/UIX compared with experiment. 
Black dashed lines show the indicated breakup thresholds for each method. The Monte 
Carlo statistical errors are shown by the light blue and yellow bands. 



Figure 5 GFMC energies using AVI 8 and AV18/IL2 compared with 
experiment. 
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Figure 6 Proton and matter densities in 4 He, 6 He, and 8 He. Analyses of 6,8 He-proton 
scattering are shown by the dashed (66) and dot-dashed (67) curves. 



Figure 7 Proton-proton distribution functions in helium isotopes. 
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Figure 8 Calculated and experimental elastic and transition form factors in 6 Li. 
Experimental data are from References 69-75. IA, Impulse approximation; MEC, 
meson-exchange contributions. 



E c . m (MeV) 

Figure 12 The total 5 , -factor for a+d radiative capture cal- 
culated for two different optical potentials with the AV18/UIX 
nuclear wave functions, compared to experimental data (91-93). 
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